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1.00 SYNPOSIS 


The basic objective of this thesis is to provide a practical method of computing the aerodynam- 
ic characteristics of slender finned vehicles such as sounding rockets, high speed bombs, and 
guided missiles. The aerodynamic characteristics considered are the normal force coefficient 
derivative, C N ; center of pressure, X; roll forcing moment coefficient derivative, O ; roll damping 
moment coefficient derivative, Q ; pitch damping moment coefficient derivative, c : ; and the drag 
coefficient, C D . Equations are determined for both subsonic and supersonic flow. No attempt is 
made to analyze the transonic region. The general configuration to which the relations are applic- 
able is a slender axisymmetric body having three or four fins. 

The fins and body are analyzed separately; and the interference between them is accounted for 
when they are combined to form the total vehicle. When necessary, both the fins and the body are 
subdivided. Both purely theoretical and semi-empirical methods are used in deriving the equations. 
Use is made of established methods when they are applicable. 

Since the high speed digital computer is a tool readily available to the contempory engineer, no 
effort has been made to linearize higher order theories or simplify complex relations. In many 
cases, the equations are derived with the computer in mind. 

The method is applied to sounding rocket vehicles and the results are compared to respective 
windtunnel and flight data. The comparison indicates that the theory predicts aerodynamic charac- 
teristics within ten percent of the experimental data. 

2.00 INTRODUCTION 
2.10 BACKGROUND 

In order to study the performance of a new or existing flight vehicle the engineer must have 
the vehicles aerodynamic characteristics. Windtunnel or flight tests can provfde the needed infor- 
mation; but they are both expensive and time consuming. Lacking either time or money the engineer 
turns to theoretical predictions. 

In this day of rockets and high speed armaments, slender finned vehicles are of considerable 
interest. A great deal of work has been done on the aerodynamics of such configurations. The list 
of contributors is exhaustive. However, a. number of works stand out and are of particular interest. 
Munk, reference 20, developed the technique of slender body theory and applied it to low speed 
airships. Tsien, reference 29, later modified Munks technique to apply it to pointed projectiles in 
supersonic flow. Jones, references 7, 8, 9, extended Munk 1 s work to finite wings and also provided 
excellent application of conical flow theory to wing lift and drag estimation. Nielsen, Pitts, Ferrari, 
and Morikawa, references 22, 3, 18 and 17, have thoroughly analyzed wing-body interference effects. 

Even though they provide valuable understanding, many of the contributions have limited appli- 
cability to practical calculations. Their inherent assumptions and approximations reduce the ac- 
curacy and reliability of their results. Those that provide accuracy are usually so complex or 
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abstract they cannot directly provide the desired aerodynamic characteristics. This thesis is an 
attempt to form concrete, useful expressions by applying the best abstract theoretical treatments 
as well as reliable empirical data. 


2.20 STATEMENT OF PROBLEM 

Determine a set of expressions for the practical calculation of the aerodynamic characteristics 
of slender, axisymmetric finned vehicles at both subsonic and supersonic speeds. The characteris- 
tics of interest are the normal force coefficient derivative, ; center of pressure, X ; roll forcing 
moment coefficient derivative, C; ; roll damping moment coefficient derivative, C p ; pitch damping 
moment coefficient derivative, C mq ; and drag coefficient, C^. 

2.30 GENERAL METHOD OF SOLUTION 

1. Divide vehicle into body and tail. 

2. Analyze the body and tail separately. Subdivide either when necessary. 

3. Analyze wing-body interference. 

4. Recombine to form total vehicle solution. 

2.40 GENERAL ASSUMPTIONS 

1. The angle-of-attack is very near zero. 

2. The flow is steady and irrotational. 

3. The vehicle is a rigid body. 

4. The nose tip is a sharp point. 

3.0 NORMAL FORCE AERODYNAMICS 

There are a number of aerodynamic coefficient derivatives that are directly dependent on the 
aerodynamic forces acting normal to the longitudinal axis of the vehicle. The normal forces acting 
on the fins are generated by four main sources; angle-of-attack; fin cant angle, roll rate, and pitch 
rate. These excitations produce the normal force coefficient derivative, C Na , at a center of pressure, 
X & y; a roll forcing moment coefficient derivative, C^ g ; a roll damping moment coefficient derivative, 
q ; and a pitch damping moment coefficient derivative, C mq respectively. The only normal force 
acting on the body is due to the angle of attack. This is the normal force at a center of pressure. 

3.10 TAIL NORMAL FORCE AERODYNAMICS 

3.11 TAIL NORMAL FORCE COEFFICIENT DERIVATIVE 

In subsonic flow, the semi-empirical method of Diederich from reference 2 gives the normal 
force coefficient derivative of a finite thin plate as; 


2 



cos r 


(3-1) 


where; 


( V 

a 1 



2 + F 


D 



C = Normal Force Coefficient Derivative of a two dimensional airfoil. 

N a0 

F = Diederich’s Planform Correlation Parameter. 

r D 

According to Diederich; 


F = * 

D 1 

— c. cos r 

277 Ao 

By the thin airfoil theory of potential flow (reference 26), 

q, ” 2 77 

"aO 

Correcting this for compressible flow according to Gothert’s rule (reference 26) 


C N 


aO 


277 

P 


Thus; 


F 


D 


/5/R 

cos r 

c 


Substituting equations (3-4) and (3-5) into (3-1) and simplifying, 


(3-2) 


(3-3) 


(3-4) 


(3-5) 



(3-6) 


This is the value for one fin. 

The value of (C^ ) for one fin in supersonic flow is computed using the method of reference 1, 
as given in Appendix A. It is essentially a strip theory which analyses two dimensional strips using 
Busemann's Third order expansion of the compressible flow equations. 


C p = K, + K, ^ + K 3 ^ + Bj rt 


(3-7) 


where; 


K, = 2 / 0 


K 


2 


(y + 1) M 4 ~ 4 /3 2 


3 



V (y + 1) M + (2T 2 - 7y - 5) M 6 + 10(7 4 1) M 4 -12 M 2 + 8 

60 > 

B, = (y 4 1) M 4 [(5 - 37) M 4 + 4 (7 - 3) M 2 + 8] J- 

J HO 

The strip values are then summed in a spanwise direction with a correction applied to those strips 
that are intersected by the fin tip mach cone. The most distinct advantage of the method in reference 
1 is that it includes a thickness distribution in the analysis. 

To apply the c. of one fin to a multi-finned vehicle the effect of fin dihedral must be con- 
sidered. Only three and four finned vehicles are analyzed, since these are the most commonly 
used. Both are analysed assuming that one fin is in a plane parallel to the flow. Reference 13 de- 
rives the effect of dihedral on the c N of two fin panels at very small angles of attack. 

a 

M, = 2 M, 00,1 A ,3 ‘ 8) 

Figure 3-1 defines the dihedral angle, a. 

A three (3) fin configuration has two effective panels at a 30° dihedral angle. Thus the tail 
normal force derivative is, 

W, = 2 ( c »„), 30 ‘ 

= 2 H 

(S.HM, « < 3 ' : 

Since the dihedral angle of the two effective panels of a four (4) fin configuration is zero, the tail 
normal force derivative is; 





(3-10) 


From equations 3-9 and 3-10 we can generalize for three and four fin vehicles. 




(3-11) 
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Figure 3-1-Fin Dihedral Angle (As Viewed From Rear) 
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Substituting 3-6 into 3-11, for subsonic flow; 



Ntt * (A f /A r ) 



(3-12) 


3.12 TAIL CENTER OF PRESSURE LOCATION 

From the thin airfoil theory of subsonic potential flow, reference 26, the center of pressure of 
a two dimensional airfoil is located at 1/4 the length of its chord from its leading edge. Thus, on 
a three dimensional fin, the center of pressure is located along the quarter chord line. By defini- 
tion, the center of pressure is also located along the mean aerodynamic chord of the fin. There- 
fore, by the above argument, the subsonic fin center of pressure is located at the intersection of 
the quarter chord line and the mean aerodynamic chord. In this treatment, the subsonic center of 
pressure is considered to remain constant with mach number. It remains to determine the length 
and position of the mean aerodynamic chord. 

The mean aerodynamic chord is defined as; 


C 


MA 




(3-13) 


Figure 3-2 shows the fin coordinate system. The generalized chord is a function of the spanwise 
coordinate. To find this function, a proportionality relation is set up. See figure 3-3. 


c 

r 

IT 



(3-14) 


It is assumed that the fin edges are all straight and form a trapizoid. From the first two terms of 
equation 3-14; 


c = c 

* r 



(3-15) 


From the first and last terms of equation 3-14; 

C L* = c L* - c S 

t r r 


Simplifying, and defining; 


K = 


(3-16) 
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Figure 3-3— Proportionality Triangle for Chord 
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we find; 


L* 


s 

l - k 


Substituting 3-17 into 3-15, 


c = c r (1 + ky) 


where, 


k = 


2_ 

v 



(3-17) 


(3-18) 


(3-19) 


Using 3-18 in 3-13; 


MA 


c 2 r s 

-t 

f Jq 


(1 + ky) 2 dy 


(3-20) 


Integrating 


c 


MA 


A, 3k 


[(1 + ks) 3 - 1 ] 


1 + k s k 2 s 2 
L 3 


(3-21) 


Substituting 3-19 into 3-21 and simplifying; 


C MA 3 


1 c ? s 


[\ 2 + \ + 1 ] 


(3-22) 


But; 

A f =i(c r+ c t )s (3-23) 

thus; 



+ c. 


[ \ 2 + K + 1 ] 
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or. 


MA 




To find the spanwise location of c A , equate 3-24 and 3-18 and solve for y. 




or 


MA 


2 2c c. 


After much algebra; 


Y - 1. 

MA 3 


c + 2c. 


C r +C t J 


The spanwise subsonic center of pressure coordinate is then; 


Yt = r * + 3 


C r + 2C «1 

. C r + C . J 


As can be seen from figure 3-1, the longitudinal coordinate is; 




MA 


But; 


l MA =y M A tanP L =__x t 


Combining equations 3-24, 3-26, 3-28, and 3-29; 






The fin center of pressure in supersonic flow is found by the method of reference 1 


(3-24) 


(3-25) 

(3-26) 

(3-27) 

(3-28) 

(3-29) 

(3-30) 

in appendix A. 
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3.13 TAIL ROLL FORCING MOMENT COEFFICIENT DERIVATIVE 

A force arises on a fin at zero angle of attack if the fin is canted at an angle with the longitu- 
dinal axis. See figure 3-4. The cant angle is an effective angle-of-attack. As can be seen in 
figure 3-5, the total force resulting from canting all the fins at the same angle is zero on three and 
four fin configurations. However, a roll forcing moment about the longitudinal axis is produced. 
Since the cant angle acts as an effective angle of attack the roll forcing moment is subsonic flow 
may be written in terms of (C^ \ ; 


^NY t (t N J 8 q A r (3-31) 

The values of (c^and Y T are found using equations 3-6 and 3-27 respectively. The rolling moment 
coefficient is defined as, 




1 

q A L 

^ r r 


(3-32) 


thus; using equation 3-26; 




N Y. 


x C N 


8 


(3-33) 


The roll forcing moment coefficient derivative is defined as; 


3Q 

Cj - — — 


(3-34) 


$ = 0 


thus; the subsonic roll forcing coefficient derivative is, 


C{ . - N M, tt 


(3-35) 


In supersonic flow, c,f is computed using the method of reference 1 as given in appendix A. 

This method determines the roll forcing moment of each strip and sums them in a spanwise direction. 


3.14 TAIL ROLL DAMPING MOMENT COEFFICIENT DERIVATIVE 

When the vehicle has a roll velocity, the superposition of the local tangential velocity of the fin 
panel on the free stream velocity produces a local angle-of-attack which is a function of spanwise 
position. The local angle-of-attack produces a local force distribution along the fin. For three 
and four fin configurations the total resulting force is zero, as shown previously. The moment pro- 
duced about the longitudinal axis opposes the roll rate (see figure 3-6) and is a roll damping moment. 
The local rolling moment at - is, 
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Figure 3-6-Roll Damping Parameters 
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dl = fF(£) 


(3-36) 


The local angle of attack at £ is, 


“o • (T) 


(3-37) 


v(0 = - P ■ 


Thus; 


a(<f) - tan 


(3-38) 


For v >>p <5 the local angle-of-attack is near zero, thus; 


“ - 


(3-39) 


In subsonic compressible flow the local normal force on an airfoil strip is found using equation 3-4. 


F ( r ) - C. q a(-) c(-) d* 
•'x 0 


(3-40) 


Where c is defined in equation 3-4. Knowing that, 

N .t o 


equation 3-18 gives 


(3-41) 


c( : r) r c 


'i-Lz± 

r L s 


(- - r ) 


... M 1 -''/ s . 

c(s)r — — (— r - 


(3-42) 
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Combining equations 3-36, 3-39, 3-40, and 3-42; 


<\i - 


f C N q 

a 0 



M 




d.f 


Simplifying; 


C N q pc r (l - \) 

dl = - a ° vg (L* + r t - C 2 d* 


(3-43) 


where L* is defined by equation 3-17. Applying equation 3-32 provides the local roll moment 
coefficient. 


dC>> 


C n Pc r Cl-M 
N a0 

K~l vs 

r r 


(L* + r { - f) ~ 2 d,~ 


(3-44) 


Integrating over the exposed fin, 


where; 


c -e 



((L* + 


••] r 2 d ' 


(3-45) 


k = - 


pc r a - M q, 

a 

A L VS 


(3-46) 


Performing the integration; 


c i = 


k s 2 

12(1 - M 


[(1 + 3 \> s 2 + 4(1 + 2 K) sr t + 6(1 + \) r 2 ] 


(3-47) 


Reintroducing k and simplifying; 


c t 


pc 


Cn 


12 BA L V 


[(1 + 3 \) s 2 + 4(1 + 2 V> sr t + 6(1 + K) r 2 ] 


(3-48) 


The roll damping moment coefficient derivative is defined as; 



Then, from 3-48, 


c r s ^ 

o = - - — [(1 + 3M s 2 4 4(1 4 2>0 sr t 4 6(1 4 K) rj] (3-50) 

p 6L 2 Ar 

T 

However, C N must be corrected for three dimensional effects using the Diederich correlation 
parameter. From the results of section 3.11, equations 3-1, 3-5, and 3-6, the subsonic roll damping 
moment coefficient derivative for N fins is, 

Nc S 

C, 1_!1 4 3\) s 2 4 4(1 + 2K) sr t 4 6(1 + K) r 2 1 (c* \ (3-51) 

p 6L 2 ' n 

r 

In order to determine the roll damping moment coefficient derivative in supersonic flow, the 
angle-of -attack distribution of equation 3-37 is introduced into the equations of reference 1 in 
appendix A. The resulting rolling moment is divided by PL r 2V to form the roll damping derivative. 
To insure that c. is truely a linear derivative, the value of n |iii must be near zero. To provide 

P 

this condition we assume 


PL 

—irk 

2V 


(3-52) 


where k is any defined constant: and, 


k << 1 


From equations 3-37 and 3-41: 


p(s 


Using equation 3-52, 


t an 


2k (s . r 


(3-53) 


For the slender, finned vehicles under consideration: 


2( s + r t ) 

L 

r 


4 


Thus 


'i ~ tan” 1 (4k) 

mn x 


(3-54) 
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To satisfy the criterion, 


a ;; tan i 


equation 3-54 dictates that, 


k 1 icr 3 

Using the equality, the angle- of-attack variation becomes; 


■002 (y * r ) 

*(y) - L ( 3 - 55 ) 

r 

This variation is introduced into the equations of reference 1 and a value of the damping C. is de- 
termined. The corresponding C; s is found by, 


For N fins; 


c t 

p 



1000 Ci 


Cj - 1 000 N Cf ( 3 - 56 ) 

p 

3.20 BODY NORMAL FORCE COEFFICIENT DERIVATIVE 

3.21 BODY NORMAL FORCE DERIVATIVE 

The axisymmetric body of a slender vehicle is generally a combination of a nose shape, circular 
cylinders and conical frustrums. The nose shape is usually conical or tangent ogival. See figure 
3-7. It is required that there be no discontinuities at the interfaces of these components. This 
allows integration over the entire length of the body. The effect of compressibility in subsonic flow 
will not be considered in this analysis. From a static stability point of view, this is a conservative 
assumption since the normal force coefficient derivative increases near mach one. 

For a slender axially symmetric body in subsonic flow, references 15 and 4 derive the steady 
state running normal load as, 


n(x) - ^V — [A(x) w(x)j ( 3 - 57 ) 

Where ACx) is the local crossectional area of the body. The rigid body downwash is a function of the 
angle-of-attack is constant in x, 


w V sin j 
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For angles-of-attack near zero 


w = Va 


Then equation 3-57 becomes, 


n (x) = pV 2 cl 


dA(x) 

dx 


From the definition of normal force coefficient derivative, 


C N (x) 


n(x) n (x) 


q A 1 , r2 . 

M r _ p V* A 

2 0 r 


Applying equation 3-60 to 3-59, 


C N (x).^^ 

A dx 


Integrating over the length of the body, 


r„ r'"dA 

v A I dx A 

r “T) r 


Thus, 


C N . [A(l 0 ) - A(0)]- 


By definition of normal force derivative, 


equation 3-63 provides, 


(3-59) 


(3-60) 


(3-61) 


(3-62) 


(3-63) 


(3-64) 


C N =L [A(l 0 ) -A(0)] (3-65) 

n A 

r 

Immediately equation 3-65 indicates that is independent of the body shape as long as the independ- 
ent of the body shape as long as the integration of equation 3-62 is valid over the body. The 
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body components being analysed are slender and have no discontinuities in their crossectional 
area or its derivatives. Thus, equation 3-65 may be applied to them. Since the nose tip is assumed 
to be a sharp point, 


A(0) = 0 

Thus, the body normal force coefficient derivative in subsonic flow is; 



(3-66) 


Reference 27 presents a very useful method for determining the normal force coefficient 
derivative of pointed bodies of revolution in supersonic flow. The second order shock expansion 
theory approximates the body by a set of tangent conical frustrums. See figure 3-8. The frustrums 
are each analyzed using an exponential series expansion for the pressure distribution 


P = - P o 


C e 


< s n s + s i 



(3-67) 


where P 0 , C, and S n are constants for a particular frustrum. By applying the boundary conditions 
at the leading edge of the frustrum and at infinity while retaining only the first series term the 
pressure distribution becomes; 


P = P - (P 

e c 


- P 2 ) c ’ 71 


(3-68) 


where; 


- = M * ~ x 2 

■' vW 2 (P c -Pj) COS Sj 
s - surface coordinate 

% 

P c = pressure on a cone having the same slope as the frustrum 
b = angle between surface and longitudinal axis 
( ) 2 = indicates the value is taken at the front edge of the current frustrum. 
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Figure 3-8- Approximate Tangent Body 
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The notation of reference 8 is used and defined to simplify correlation between this presentation 
and that reference. The value of (ap/?ts) 2 is determined by an approximate solution of the general- 
ized characteristic equation in axially symmetric flow, 

"P db_ _ 1 3P - K ( as> + S i n /i sin A (3-69) 

as Js cos fi SCj cos ix \ aCj r J 


where 

>. = 2 i£ 

sin 

C = First Characteristic Coordinate 
= Mach Angle 
r = local radius 
y = ratio of specific heats 

The solution is; 



(3-70) 


where; 

y .M 2 

B = —1 

2 (M 2 - n 


1 1 

i * 1 

f '- v 


^ 2 j 

M 


y + l 



2 


( ) indicates the value is taken at the trailing edge of the preceeding frustrum. 
/ \ indicates the value is taken at the trailing edge of the current frustrum. 

k 3 

The value of (dP/.-s) is just ( hp ds). of the previous frustrum. 

i 3 



(3-71) 


By applying equations 3-68, 3-70, and 3-71 to each element in succession, the pressure distribution 
over the tangent body is determined. The local nondimensional normal load is given by; 
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2 


d(P/P.) 


cos p dp 


(3-72) 


A - n( _ x) 
a q d 


y M, 


“I" 

i 77 *'0 


da 


where 

<p = Circumferential Angle, 
d = Local diameter. 

The derivative of the pressure distribution is evaluated by differentiating equation 3-68 with respect 
to angle of attack; and then applying physical restraints to simplify the result. 


dCP/PJ 

da 


a - e~ v ) 


d ( p e /p *> 

da + 0 A, 


d(P, /P,) 

da 


Using 3-73 in 3-72 and simplifying, 


A = (1 - e" 77 ) (tan S') 




(3-73) 


(3-74) 


Where "tcx” indicates the value of an equivalent cone. The running load distribution is determined 
by successive application of equation 3-74 to the tangent body elements. The value of Aj, for the 
initial cone is just, 

A - (tan & v ) (c \ (3-75) 

'“/lev 

Where "v" and "tcv” are the values at the initial equivalent cone vertex. The normal force coefficient 
derivative is determined by the longitudinal integration of the running load distribution. 



(3-76) 


where ; 

A b = Body Base Area. 

1 = Body Length. 

The initial cone and equivalent cone characteristics are obtained from references 24 and 25; and are 
presented in tables 3-1 thru 3-3. 
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. Table 3-1-Values of M at Initial Cone Surface 



2.5 

5.0 

7.5 

10.0 

12.5 

15.0 

17.5 

20.0 

22.5 

25.0 

27.5 

30.0 

1.5 

1.4866884 

1.4579318 

1.4197316 

1.3748579 

1.3248992 

1.2707392 

1.2127162 

1.1506488 

1.0837313 

1.0100703 



1.75 

1.7340335 

1.7004643 

1.6568043 

1.6064462 

1.5513276 

1.4925484 

1.4306418 

1.3657361 

1.2976505 

1.2259024 

1.1495789 

1.0669388 

2.0 

1.9808615 

1.9415254 

1.8912335 

1.8340453 

1.7721933 

1.7068837 

1.6386764 

1.5677572 

1.4940884 

1.4174680 

1.3375202 

1.2535953 

2.5 

2.4729472 

2.4193976 

2.3528628 

2.2788768 

2.2001554 

2.1179422 

2.0327606 

1.9448503 

1.8543543 

1.7613748 

1.6659616 

1.5680718 

3.0 

2.9627600 

2.8913832 

2.8048038 

2.7101402 

2.6103929 

2.5067511 

2.3997895 

2.2899646 

2.1777409 

2.0635830 

1.9479054 

1.8310067 

3.5 

3.4500939 

3.3571597 

3.2467499 

3.1275136 

3.0025569 

2.8731152 

2.7400668 

2.6043443 

2.4669130 

2.3286833 

2.1904152 

2.0526629 

4.0 

3.9347700 

3.8164543 

' 3.6783916 

3.5305680 

3.3761469 

3.2166972 

3.0537628 

2.8889908 

2.7239783 

2.5600783 

2.3983186 

2.2393739 

4.5 

4.4166246 

4.2690580 

4.0994611 

3.9189007 

3.7307750 

3.5374323 

3.3414321 

3.1453333 

2.9513677 

2.7612386 

2.5760876 

2.3965344 

5.0 

4.8955305 

4.7148078 

4.5097079 

4.2921847 

4.0662668 

3.8355932 ! 

3.6040616 

3.3752078 

3.1517935 

2.9356870 

2.7279347 

2.5289133 

6.0 

5.8440686 

5.5852474 

5.2968787 

4.9927832 

4.6801702 

4.3668118 

4.0597311 

3.7640785 

3.4829712 

3.2178167 

2.9687934 

2.7352580 

7.0 

6.7796857 

6.4269172 

6.0384888 

5.6318010 

5.2204850 

4.8178049 

4.4335119 : 

4.0731395 

3.7388047 

3.4303400 

3.1462500 

2.8843548 

8.0 

7.7018785 

7.2390804 

6.7337145 

6.2103924 

5.6922672 

5.1981506 

4.7389543 ! 

4.3187765 

3.9373126 

3.5918831 ! 

3.2787858 

2.9940481 

10.0 

9.5046400 

8.7725147 

7.9851813 

7.1977580 

6.4568340 

5.7867581 

5.1933727 

4.6723717 

4.2153623 

3.8131019 

3.4569163 

3.1391940 

12.0 

11.250395 

10.182053 

9.0582962 

7.9858490 

7.0291308 

6.2043028 

5.5019835 

4.9042678 

4.3926921 

3.9510537 

3.5659952 

3.2267626 

15.0 

13.759139 

12.063489 

10.368285 

8.8727901 

7.6324201 

6.6228697 

5.7997803 

5.1216448 

4.5552592 

4.0753516 

3.6629457 

3.3037471 

20.0 

17.637850 

14.613515 

11.919518 

9.8177450 

8.2281047 

7.0145241 

6.0679293 

5.3119812 

4.6946790 

4.1802823 

3.7437969 

3.3673365 




Table 3-2-Values of — - 

\ da. 


a = 0 


for Equivalent Cones 


\ s 
M \ 

2.5 

5.0 

7.5 

10.0 

12.5 

15.0 

17.5 

20.0 

22.5 

25.0 

27.5 

30.0 

1.5 

1.9770483 

1.9336127 

1.8821566 

1.8280992 

1.7735619 

1.7187816 

1.6630926 

1.6057182 

1.5462941 

1.4852936 



1.75 

1.9731742 

1.9258447 

1.8737533 

1.8224578 

1.7730646 

1.7243971 

1.6746400 

1.6223425 

1.5666956 

1.5075064 

1.4451701 

1.3808820 

2.0 

1.9690287 

1.9180672 

1.8661051 

1.8183423 

1.7742540 

1.7310117 

1.6859216 

1.6372641 

1.5842135 

1.5265736 

1.4645815 

1.3988695 

2.5 

1.9602967 

1.9036740 

1.8548454 

1.8161039 

1.7820973 

1.7472553 

1.7083329 

1.6638860 

1.6134668 

1.5571606 

1.4953496 

1.4286279 

3.0 

1.9514206 

1.8917127 

1.8493550 

1.8202950 

1.7941549 

1.7645092 

1.7287097 

1.6859461 

1.6362061 

1.5798551 

1.5174588 

1.4497109 

3.5 

1.9426977 

1.8826075 

1.8487600 

1.8282076 

1.8074148 

1.7806466 

1.7461812 

1.7037859 

1.6537994 

1.5968132 

1.5335317 

1.4647351 

4.0 

1.9343213 

1.8763907 

1.8516843 

1.8377630 

1.8203259 

1.7949749 

1.7608321 

1.7181330 

1.6674853 

1.6096523 

1.5454397 

1.4756774 

4.5 

1.9264473 

1.8728601 

1.8568202 

1.8477350 

1.8322685 

1.8074272 

1.7730368 

1.7296967 

1.6782310 

1.6195188 

1.5544309 

1.4838272 

5.0 

1.9191785 

1.8715909 

1.8632046 

1.8574765 

1.8430453 

1.8181556 

1.7832011 

l '.7390731 

1.6867491 

1.6272149 

1.5613461 

1.4900257 

6.0 

1.9066985 

1.8741557 

1.8772933 

1.8752450 

1.8611884 

1.8352940 

1.7988073 

1.7530357 

1.6991506 

1.6381839 

1.5710540 

1.4986224 

7.0 

1.8972025 

1.8808308 

1.8911403 

1.8902959 

1.8754007 

1.8479983 

1.8098978 

1.7626404 

1.7074620 

1.6453960 

1.5773421 

1.5041237 

8.0 

1.8906735 

1.8893943 

1.9037986 

1.9027991 

1.8865235 

1.8575153 

1.8179511 

1.7694507 

1.7132534 

1.6503536 

1.5816157 

1.5078364 

10.0 

1.8854155 

1.9075860 

1.9248544 

1.9216225 

1.9022129 

1.8703588 

1.8284727 

1.7781460 

1.7205191 

1.6564935 

1.5868584 

1.5123524 

12.0 

1.8876297 

1.9243231 

1.9408134 

1.9345246 

1.9122898 

1.8782443 

1.8347389 

1.7832191 

1.7246910 

1.6599792 

1.5898124 

1.5148790 

15.0 

1.8980460 

1.9449893 

1.9576949 

1.9470431 

1.9215589 

1.8852675 

1.8402027 

1.7875704 

1.7282369 

1.6629184 

1.5922831 

1.5169846 

20.0 

1.9201707 

1.9690214 

1.9743436 

1.9583547 

1.9295376 

1.8911343 

1.8446777 

1.7910947 

1.7310758 

1.6652577 

1.5942431 

1.5186494 




Table 3-3-Values of Surface Pressure Coefficient for Equivalent Cones 



2.5 

5.0 

7.5 

10,0 

12.5 

15.0 

17.5 

20.0 

22.5 

25.0 

27.5 

30.0 

1.5 

.0123382 

.0396615 

.0773864 

.1237983 

.1781167 

.2400313 

.3095450 

.3869791 

.4731538 

.5699888 



1.75 

.0114414 

.0363287 

.0704289 

.1122701 

.1612014 

.2169225 

.2792900 

.3482459 

.4238192 

.50619022 

.59585744 

.69405243 

2.0 

.0107862 

.0339457 

.0655812 

.1044583 

.1500854 

.2022339 

.2607607 

.3255312 

.39641336 

.47330179 

.55616008 

.64511954 

2.5 

.00982981 

.03058576 

.05895952 

.09413081 

.1358608 

.1840343 

.23852578 

.29915906 

.36571047 

.43791904 

.51550425 

.59818912 

3.0 

.009 1 3 5 41 

.02824024 

.05451919 

.08747518 

.1270239 

.17310135 

.22558051 

.28426241 

.34888446 

.41912862 

.49463659 

.57503078 

3.5 

.00859242 

.02647476 

.05130225 

.08281755 

.12101992 

.16585924 

.21719882 

.27481810 

.33842506 

.40766937 

.48215440 

.56145326 

4.0 

.00814929 

.02508810 

.04886536 

.07939344 

.11671212 

.16077021 

.21141454 

.26840424 

.33142494 

.40010296 

.47401780 

.55271448 

4.5 

.00777855 

.02396847 

.04696202 

.07678924 

.11350386 

.15704447 

.20724252 

.26383657 

.32649601 

.39*482976 

.46839852 

.54673311 

5.0 

.00746199 

.02304657 

.04544267 

.07475691 

.11104466 

.15423167 

.20413070 

.26046646 

.32289183 

.39100288 

.46435008 

.54245194 

6.0 

.00694671 

.02162337 

.04318936 

.07182851 

.10757999 

.15033841 

.19988673 

.25592332 

.31807966 

.38593639 

.45902778 

.53685790 

7.0 

.00654427 

.02058479 

.04162143 

.06985826 

.10530922 

.14783838 

.19720408 

.25308627 

.31510472 

.38282933 

.45578647 

.53347449 

8.0 

.00622092 

.01980211 

.04048580 

.06847199 

.10374427 

.14614133 

.19540407 

.25120009 

.31314047 

.38078911 

.45366854 

.53127090 

10.0 

.00573592 

.01871810 

.03898810 

.06670277 

.10179238 

.14406008 

.19322118 

.24893118 

.31079304 

.37836321 

.45115963 

.52867228 

12.0 

.00539243 

.01802077 

.03807635 

.06566456 

.10067415 

.14288515 

.19200207 

.24767463 

.30949966 

.37703244 

.44978975 

.52725561 

15.0 

.00503756 

.01736467 

.03726273 

.06476571 

.09972302 

.14189705 

.19098631 

.24663061 

.30843098 

.37593703 

.44866381 

.52609497 

20.0 

.00468137 

.01677876 

.03657736 

.06403097 

.09895953 

.14111270 

.19018534 

.24581285 

.30759370 

.37507865 

.44778536 

.52519347 




3.22 BODY CENTER OF PRESSURE LOCATION 


As in the subsonic body normal force coefficient derivative analysis, it is required that there 
be no discontinuities between the body component shapes in order to allow integration over the en- 
tire body. 

By definition, the pitching moment of the local aerodynamic force about the front of the body 
(x = 0) is; 


% 

m(x) 


xn (x) 


(3-77) 


The pitch moment coefficient is defined as; 


C 

m 


m 

EjT~L 

^ r r 


(3-78) 


thus, 


Since, 


xn(x) 
q A L 

M r r 


Cn 



equation (3-79) becomes 


C 

m 


(X) 


x Cj^Cx) 

L 

r 


(3-79) 


(3-80) 


(3-81) 


The value of the subsonic local normal force coefficient is given in equation 3-61. Using it in 
equation 3-81. 


c (x) 


2 ax dA(x) 

A L dx 


(3-82) 


Integrating over the length of the body, 


C 



(3-83) 
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This can be integrated by parts; 



dx 


giving; 


C 

m 


2 a 

rr 

r r 


f l ° 

-ig A C 1 o - A(x) dx 


(3-84) 


But, by definition, the integral term is the body volume, therefore; 


C 

m 


2a 

A L 

r r 


-v 


By definition, the pitch moment coefficient derivative is; 

3C I 

n m i 


(3-85) 


(3-86) 


Then; from equation 3-85, 



2 

A L 

r r 


n 0 A (V - \ ] 


The center of pressure is determined by, 



Using equations 3-65 and 3-87 in 3-88 and simplifying yields, 


(3-87) 


(3-88) 


x - l 


(3-89) 
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The volume formulae for cones, conical frustrums and circular cylinders are well known and 
do not merit repetition in this paper. However, the volume formula for a tangent ogive is not 
readily available and is given herein for the readers convenience. The volume of a tangent ogive 
nose shape is determined in reference 12 as; 



In supersonic flow, the method of reference 27 determines the pitch moment coefficient derivative as 


C 



A rx dx 


(3-91) 


Where \(x) is found using equation 3-74. The center of pressure location of the body in supersonic 
flow is then determined by the use of equations 3-76 and 3-91 in equation 3-88, 


X 


B " 



(3-92) 


3.30 TAIL-BODY INTERFERENCE EFFECTS 

The effect of interference flow between the tail and the body on the normal force aerodynamics 
of a slender, finned vehicle can be expressed in terms of interference factors (reference 23). 


(u., ■ to K "» 

(3-93) 

to<T> = to K *' , > 

(3-94) 

(to,.; to 

(3-95) 


The interference factors are determined by applying slender body theory to the flow over the total 
vehicle. 
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3.31 TAIL IN THE PRESENCE OF THE BODY 

Because the fins under consideration are of low aspect ratio, the value of K T(B) is available 
from references 7 and 17. 


K 


T( B ) 






(3-96) 


Reference 23 determines that the longitudinal tail center-of-pressure is independent of the inter- 
ference flow. 


X T(B> = X T 


(3-97) 


3.32 BODY IN THE PRESENCE OF THE TAIL 

The only portion of the body affected by interference flow is the portion between and behind the 
fins. Subsonically, slender body theory, reference 3, determines; 


K 


B(T) 



(3-98) 


where k t(B) is given by equation 3-96. In supersonic flow, the interference effect on the body can 
be strongly affected by the fin tip mach cone. The fin tip mach cone intersects the body when, 

/3/R (1 + V) + 1 j >4 (3-99) 


Until this condition is reached equation 3-98 is still valid in supersonic flow. Once the fin tip 
mach cone impinges on the body, the function becomes quite complex and is more easily presented 
in the form of curves. Figure 3-9 shows the variation of the corrected interference factor. Using 
the result of this curve, K^ (T) ; 
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(3-100) 



K' 

B(T) 


(i • *) (r - n 


By the extension of the fin lifting line into the body, reference 23 derives the subsonic variation 
of the center-of-pressure location of the tail interference flow over the body. 


X - 1 - _L, 

b(t> T 4 


cosh 1 " 


Sr r 




cosh" 


‘(tth-i 


tan fj /4 (3-101) 


Where, s > r and /3AR - 4. In addition, slender body theory is applied subsonically to find X B at 
;*AR= 0. Slender body \ (T) values are obtained at intermediate sfR values only for \ = 1. For other 
values of \ and sweep the 7 B(T) curves are interpolated between the lifting line values and the 
slender body value at 5AR= 0. The shape of the interpolated curves is assumed similar to the K = l 
curve. Figure 3-10 shows the full subsonic *X B(T) variation. In supersonic flow, the center-of- 
pressure of the interference flow over the body is found by applying equation 3-92; 


X 


B(T) 



d + 1 


T 


(3-102) 


Reference 23 gives the body interference pitch moment derivative as derived from slender body 
theory; 



8 (m/3) 372 

+ 1) 



X cos 



(3-103) 


This is for fins with supersonic leading edges, m3 > 1. For subsonic leading edges, m3 < 1 and, 



(3-104) 
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K B{ T )- K B[T + 1 ^ ^ 



<c>v 


Figure 3-9 — Supersoni c Kg^yj 
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(d)NO MIDCHORD SWEEP, X 


EFFECTIVE ASPECT RATIO, /3 /R 
=0. (e)NO MIDCHORDSWEEP, X= y- (f )NO MIDCHORD SWEEP, X = l 


Figure 3-10 (Continued) — Subsonic X^y^ 






B(T) 






Where the coordinate system has been given in figure 3-2. Xg T) is determined by using equation 
3-100 along with either 3-103 or 3-104 inequation 3-102. The result is shown in figure 3-11. 

3.33 CANTED TAIL IN THE PRESENCE OF THE BODY 

The effect of fin cant at zero angle-of-attack is determined in a manner similar to the effect 
at an angle-of-attack alone. Reference 23 determines the corresponding interference factor, 



(3-105) 


No center of pressure value need be considered for the canted fin in the presence of the body 
since only the roll moment derivative is considered. 


3.40 TOTAL NORMAL FORCE AERODYNAMICS 

The total normal force coefficient derivative is the sum of all the components. 

c n --M ♦ M 

\ -fB \ a Ar ( B ) \ VB(T) 

The overall center of pressure is the solution of the moment sum of the components. 

^ _ Xb Kl + Xt(B) (Sl ( B) + Xb(T) Mg, 


(3-106) 


(3-107) 


The value of total Y is zero. The total Q value is determined from equation 3-95. 
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3.41 PITCH DAMPING MOMENT COEFFICIENT DERIVATIVE 

When the vehicle has an angular pitch velocity, the superposition of the local pitch velocity at 
the fins on the free stream velocity produces a local angle of attack which is constant across the 
fins. See figure 3-12. The moment produced by the induced force is; 
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V ANGLE -OF- ATTACK 

Figure 3-12-Pitch Domping Parameters 
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m F(Ax) 


(3-108) 


where; 


But, 


and, 


T(B) 


F = fc N \ A r q (3-109) 

\ J /T(B) 

* 

a' = tan"' (J) 


For V » v, 


(3-110) 


the local pitch velocity is, 


V = q(Ax) 


(3-111) 


Combining equation 3-108, 3-109, 3-110, and 3-111 and simplifying; 


m 



(AX ) 2 A, 


(?) 


<1 77 


(3-112) 


Knowing the definition of pitch moment coefficient, equation 3-78, 


C 

m 




(3-113) 


The pitch damping moment coefficient derivative is defined as, 


C 

mq 



(3-114) 


Thus, using equation 3-113 in 3-114. 
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(3-115) 


Equation 3-123 is equally applicable to three and four fin configurations since the tail normal force 
derivative is corrected for dihedral effects. 


3.42 ROLL DAMPING IN THE PRESENCE OF THE BODY 

Any tail body interference factor may be expressed as the ratio, 


K 



a.(') c(~) dc 


ol 0 (£) c(') d ' 


(3-116) 


(See figure 6) where; 

x ('■) = the angle -of- attack variation over the fin in the presence of the body 
x (;) = the angle-of-attack variation over the fin alone. 

Reference 23 gives the value of a t (£) for a tail-body combination at an angle of attack, a, as; 




a 



(3-117) 


where : B is the angle of attack at the body surface. For the case of the tail-body combination at an 
angle of attack, 


x B = a = a 0 (O 


However, for a rolling vehicle at a zero angle of attack, equation 3-39 gives; 




P£ 

v 


and 



Then, equations 3-39, 3-117, and 3-118 give, 



(3-118) 


(3-119) 
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Therefore, the roll damping interference factor, k 


, becomes 


R(B) 


Simplifying; 



Integrating equation 3-121 using the chord variation given by equation 3-42 yields, 


(3-120) 


(3-121) 


k R(B) * + 


1 


1 n 


(r + 1 ) (r - K) (1 - >) ( r 3 - 1 ) 

' “ 2 “ 3(t - 1) 


( 3 - 122 ) 


The roll damping moment coefficient derivative can then be expressed as, 



( C ^p) k R(B) 


(3-123) 


Where k is given by equation 3-122. 

R( B) 


4.00 AXIAL FORCE AERODYNAMICS 

The only coefficient dependent on the aerodynamic forces acting parallel to the longitudinal 
axis of the vehicle is the drag coefficient. Drag forces arise from two basic mechanisms, skin 
friction and pressure distribution. Since both mechanisms, and their interactions, are extremely 
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complex, most practical theoretical drag computations must be tempered with statistical and 
empirical data. The majority of the equations derived in this section are either theoretical gen- 
eralizations of statistical-empirical data or direct curve fits of experimental results. 


4.10 SKIN FRICTION DRAG 

The skin friction coefficient is defined as the skin friction drag of a plane surface divided by 
the dynamic pressure and the wetted area of the surface. 


4.11 INCOMPRESSIBLE FLOW SKIN FRICTION COEFFICIENT 

The skin friction coefficient for incompressible flow can be expressed in terms of the reynolds 
number, see figure 4-1. For laminar flow, their relation has been established both experimentally 
and by theoretical boundary layer considerations, reference 6. 


1.328 

f " /F 

e 


(4-2) 


For reynolds numbers greater than about 5 x 10 s the boundary layer flow becomes transitional and 
then turbulent. No good theoretical solutions have been found for turbulent flow. Figure 4-1 shows 
the experimental data for turbulent skin friction. Schoenherr applied the theoretical boundary layer 
work of von Karman to the experimental data and established the turbulent flow relationship; 


log (R, C f ) 


0.242 



(4-3) 


Hama formed a simpler formula which agrees with equation 4-3 within ±2%. 

C f = 1 (4-4) 

(3.46 log R e - 5.6) 2 

The form of the transitional curve may be approximated by subtracting an increment from equation 
4.4. The increment has the form 


AC f = h- (4-5) 

e 

Based on experimental data Prandtl indicated thatk = 1700. This fact, as well as other possible 
transitional curves are shown on figure 4-1. Combining equations 4-4 and 4-5, 
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Figure 4-1 —Skin Friction Coefficient 




Table 4-1 

Approximate Surface Roughness Heights of 
Physical Surfaces 


Type of Surface 

Approximate 

microns 

R s in mils 

Surfaces like that of a "mirror" 

0 

0 

Surface of average glass 

0.1 

0.004 

Finished and polished surfaces 

0.5 

0.02 

Aircraft-type sheet-metal surfaces 

2 

0.1 

Optimum paint- sprayed surfaces 

5 

0.2 

Planed wooden boards 

15 

0.6 

Paint in aircraft-mass production 

20 

1 

Steel plating - bare 

50 

2 

Smooth cement surface 

50 

2 

Surface with asphalt-type coating 

100 

4 

Dip- galvanized metal surface 

150 

6 

Incorrectly sprayed aircraft paint 

200 

8 

Natural surface of cast iron 

250 

10 

Raw wooden boards 

500 

20 

Average concrete surface 

1000 

40 
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1 


1700 


(4-6) 


C r = 


(3.46 log R - 5.6) 2 


for R > 5 * 10 5 . 

e 

The above considerations are for surface roughness sizes that are submerged in the laminar 
sublayer. For appreciable rough surfaces the turbulent skin friction coefficient is shown in figure 
4-2. The "smooth" curve is followed until the critical reynolds number corresponding to the rough- 
ness is reached. 


(R )„ 



(4-7) 


For reynolds numbers higher than the critical value, the skin friction coefficient can be considered 
independent of reynolds number, reference 6. 


C 


f 



(4-8) 


The flow is considered to be turbulent throughout the reynolds number range since the roughness 
trips transition even at low reynolds numbers. Table 4-1 is a list from reference 6 of the approxi- 
mate average roughness heights for physical surfaces. 


4.12 COMPRESSIBLE SUBSONIC VARIATION 

In compressible subsonic flow reference 6 indicates that the laminar skin friction coefficient 
remains unchanged while the "smooth" turbulent value varies as, 


C f = C f (1 - .09 M 2 ) 

c 

and the change for turbulent flow with roughness is, 

C f = C f (1 - .12 M 2 ) 


(4-9) 


(4-10) 


4.13 SUPERSONIC VARIATION 

Theoretical solutions in reference 16 indicate that the supersonic variation of the Laminar skin 
friction coefficient is; 


C - 1 (4-11) 

f c (1 + .045 M 2 ) 1/4 

This has been confirmed by experimental data, reference 6. Experimental data, reference 6, for 
the variation of "smooth" turbulent C f in supersonic flow indicates a strong dependence on the wall 

c 
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heat transfer as well as mach number, figure 4-3. The wall temperature alternatives shown are 
the two possible extremes. The actual condition would probably be closest to the recovery tempera- 
ture. The curve fit of figure 4-3 is, 


(1 + k M 2 )- 58 


(4-12) 


where k = .15 for no heat transfer and k = .0512 for cooled wall. 


According to reference 10, the turbulent sktn friction with roughness drops off according to the 
relation, 




1 + .18 M 2 


(4-13) 


However, this equation must be used with care. It is apparent the value of C f calculated for the 
flow with roughness should never be less than the corresponding smooth value. 


4.14 DRAG COEFFICIENT 

Once the skin friction coefficient is determined, the skin friction drag coefficient, as based on 
reference area, is found by, 


C 


D 




(4-14) 


For the fins, the wetted area of each fin is twice it r s planform area. Thus, the tail skin friction 
drag coefficient is, 


C D f =2NC fc 



(4-15) 


The body skin friction coefficient must be corrected for the fact that the body is not a flat plate, 
reference 6. 



(4-16) 


4.20 PRESSURE DRAG 

4.21 TAIL PRESSURE DRAG 

The pressure drag on a fin in subsonic flow arises from three basic contributions: a finite 
radius leading edge; thick trailing edge; and the overall fin thickness. 
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Figure 4-3— Supersonic Cj Variation 
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A finite radius leading edge can be considered a circular cylinder in cross flow that effectively 
has no base drag. Figure 4-4 from reference 6 shows the variation of the total drag coefficient of 
an infinite circular cylinder perpendicular to the flow and the variation of its base drag. The dif- 
ference between these two curves is the drag coefficient of the leading edge. This difference is 
plotted in figure 4-5. 

Subsonically up to M = .9 this curve is fitted by, 

AC d = (1 -M 2 )~- 417 - 1 (M < .9) (4-17) 

Between .9 < M < 1, the curve is essentially linear. 

ACjj = 1-1.5 (M - .9) (.9<M<1) (4-18) 

Supersonically, the curve may be approximated by a high order polynomial, reference 6. 


£C p = 1.214 - l?- 02 . + (M > 1) (4-19) 

M 2 M 4 M 6 

From cross flow theory, reference 6, the drag coefficient of a cylinder inclined to the flow is, 


~D/ 


s 3 r 


(4-20) 


The frontal area of a typical leading edge is; 


\e - ^ l le r L 


the length can be expressed in terms of the span, 


thus; 


Lle " cos r. 


A LE 2Sr L 


A A cos r. 

r r L 


From equations 4-20 and 4-21, the total leading edge drag of N fins becomes, 

\ = 2N m c ° s2 ^ ^ 


(4-21) 


(4-22) 


where AC D is obtained using equations 4-17 thru 4-19. 

The trailing edge drag is two dimensional base drag and can be expressed in terms of the cor- 
rected skin friction coefficient, reference 5, for incompressible flow. 


0.135 


% 3 




(4-23) 


50 



CYLINDERS BETWEEN WALLS 

* ARC, TUNNEL TESTS 

0 AVA, WIRES i TUNNEL 

» NACATUNNEL 

* NACA, FLIGHT TEST 

* BRITISH, STANTON 

* DVL, SUBSONIC 

▼ JUNKERS, TUNNEL 

* NACA SUBSONIC 

BASE DRAG 
/' SUBSONIC CYLINOER 
n SUBSONIC CYLINDER 

* TRANSONIC CYLINDER 

* SUBSONIC SUPERSONIC 

1 SPHERE (CORRECTED) 

> NAVORD, SPHERE 


O.l 0.2 0.3 0.4 0.6 0.8 I 2 3 4 6 8 10 

MACH NUMBER 



Figure 4-4-Cylinder Drag Variation 
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where; 



(4-24) 


For subsonic compressible flow, this may be corrected using the prandtl factor for swept wings, 
reference 6, 


P 


r 


l 

y/ 1 - M 2 cos 2 T c 


(4-25) 


However, in order to maintain a finite value at 

M 2 cos 2 P c - 1 

the prandtl factor is corrected to, 


P - r 1 (4-26) 

y/^K - M 2 cos 2 

where K is a constant to be determined from the value of C D at M = 1. Also, for compressible flow, 

b t 

equation 4-24 becomes, 


C 


2 C r 



(4-27) 


After correcting to the reference area, the subsonic trailing edge drag coefficient for N fins is, 



(4-28) 


Reference 11 gives a curve for the two dimensional base pressure coefficient at supersonic speeds 
figure 4-6. This is essentially the base drag coefficient refered to the base area. An analytic ex- 
pression which remains well within the cross hatching is, 

„ 1 - .52 IT 119 


But, this must be corrected for possible boattailing of the fin, reference 6. 


p - 

1 - .52 M-> 19 


*B - 

/t \ 2_ 



1 + 18C f Ur) 

M 2 


c \ h ,/ 



(4-29) 
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For N fins, the supersonic trailing edge drag coefficient, as corrected to the reference area is then, 



NCI - -52 M _1 

..9) | 

\ A r / 


00 

n 

nV 

M 2 


L 

\uj 



In order to determine the constant K, equations 4-28 and 4-30 are set equal at M = 1. 


(4-30) 



where the superscript * indicates sonic values. Solving forK, 



(4-31) 


This is used in equation 4-28 to evaluate the subsonic compressible trailing edge drag. 


This method of maintaining finite values of the prandtl factor at mach one is used throughout 
the remaining derivations. 

From reference 6, the drag due to thickness of two dimensional airfoil in compressible sub- 
sonic flow can be expressed as, 


C„ C. 
d tt f c 


4 (I 


120 


(;)' 


(Vl - M 2 )' 


For a three dimensional fin, based on the reference area this becomes, 


CL - C, — 


f c \A 


120 


4|—| f 
c 


ft! 


^ 1 - M 2 cos 2 I 


(4-32) 


(4-33) 
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Correcting it for a fin sweep and introducing the correction constant in the prandtl factor, 






( K - M 2 cos 2 ) 3 2 


(4-34) 


The value of the thickness drag of finite fins at mach one has been determined and presented in 
reference 6. A curve fit of the result is, 


C* rl.ls/'-ll) [l.61 + .f- yVf - 1.43) 2 + .578] 

TT VC / L J 

where 



Solving 4-34 for K, 



(4-35) 


(4-36) 


Equation 4-35 is used in 4-36 to determine K. The result is used in the following equation to de- 
termine the subsonic tail thickness drag. For N fins, equation 4-34 is; 


4N C 




30 | 

rt 'V 2 r ~ 1 

COS* 1 

K C r) 

. _ . 

(A 

- M 2 cos 2 r t j | 


(4-37) 


At supersonic speeds the effect of a finite radius leading edge is found by using equation 4-19 in 

equation 4-22. The trailing edge drag is calculated using equation 4-30. The thickness or wave 

drag, C D , is computed by the method of reference 1 in appendix A. 
w T 

The overall drag on the tail is the sum of the different contributions. In subsonic flow 




(4-38) 
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At supersonic speeds, 



(4-39) 


Equation 4-38 uses the results from equations 4-15, 4-22, 4-28, and 4-37 respectively. Equation 
4-39 uses equations 4-15, 4-22, 4-30, and results of appendix A respectively. 


4.22 BODY PRESSURE DRAG 

The body pressure drag may be broken down into two contributions, for body and base. Accord- 
ing to reference 6, the forbody drag in subsonic compressible flow may be written, 

6 A w B c f 

d r 

c - 

D P /, V ' 2 

f3 A r (/iTtf) 

Again applying the correction constant to the prandtl factor, 


D p f| A r (K - M 2 y 6 

From experimental data, reference 19, the forbody drag values of both ogives and cones are known, 
figure 4-7. In the range of values that are of interest, 1 1 f B i 10, these curves may be fitted by 
the following equations. 


(4-40) 


(4-41) 





.88 

+ -7) 1 29 


(4-42) 



.88 


(f N t n 2 


Solving equation 4-41 for k; 


K = 1 


6A W C* f 

*B c 

f a \ c ;~ 


5/3 


(4-43) 


(4-44) 
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Then 


K = 1 




6.82 A w C* f (f„ + .7)« 


f B \ 


(4-45) 


K 

ogive 


1 * 


‘6.82 a w c: (f t n 2 - 22 

B f c N 


*1 \ 


(4-46) 


The subsonic compressible forbody pressure drag is found by using equations 4-45 or 4-46 in equa- 
tion 4-41. For the nose and f should be used. For a conical shoulder A w and f s are 

W N a n a s a a 

proper. 

The supersonic value is determined by the second-order shock-expansion method, reference 
27, as explained in section 3.21 of this thesis. The axial component of the pressure distribution 
is integrated around each element and then over the length of the body. 


n ’T 

Psin cl.-Hx (4-47) 

. , 


where 

P comes from equation 3-69 

o = the local angle between the surface and the longitudinal axis 
,/) = Circumferential Angle. 

The Base drag of a body of revolution in incompressible flow is presented in reference 5 as, 


where, 


0.029 





Correcting this to the reference area and for subsonic compressibility. 




(1 - M 2 ) 


(4-48) 


(4-49) 
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Introducing the correction constant into the prandtl factor term. 


Cd 


B 



(4-50) 


From reference 11, for no fins near the body base, or for fins of zero thickness, the body base 
pressure coefficient at mach one is 


P* r 0.185 


At mach one with h f /c p = 0.1, 


Pg = 0.30 


See figure 4-8. 

The results of reference 11 indicate that P B is roughly a linear function of fin trailing edge 
thickness at a constant mach number. Therefore, 


Solving equation 4-50 for K, 



(4-51) 


(4-52) 


Using equation 4-51 in 4-52 and simplifying 



(4-53) 


The subsonic compressible body base drag is computed by using equation 4-53 in 4-50. 

References 11 and 6 both present curves of the base drag in supersonic flow. Figure 4-8 is 
from reference 11. The free flight test curve is most representative of the variation. It is an ex- 
ponential over an initial range and an inverse square over the remainder of the mach numbers. To 
find the exponential, the free flight curve in figure 4-8 may be fitted in terms of the value at mach 
one. 
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(4-54) 


C C* [0.88 t 12 e~ 3 ' 58(M 1 } 1 (1 ' M < M rr ) 

d b d b 

where M is the point where the two variations cross. Also, reference 5 gives the inverse square 
variation as, 

.7 


If it is assumed the ratio of the value of at M = M cr to C* remains constant for all curves, (inspec- 
tion of the curves in reference 6 indicates this is a good assumption) then the value of the critical 
mach number can be found, 



0.892 (4-56) 



The value of C* is given by equation 4-51. In either subsonic or supersonic flow, 



For subsonic values, equation 4-41 and 4-50 are used. In supersonic flow equations 4-47 and 4-54 
or 55 are used. 


4.23 TOTAL DRAG 

The drag of the total vehicle is the sum of the tail and body contributions. 


No interference drag effects are considered. 

5.00 COMPARISON WITH EXPERIMENTAL DATA 

In an attempt to obtain as broad a range as possible without being voluminous, each characteris- 
tic is compaired to the data from a different configuration. In some cases, the choice of configuration 
was dictated by the availability of the appropriate data. 

5.10 NORMAL FORCE COEFFICIENT DERIVATIVE AND CENTER OF PRESSURE 

The windtunnel data of the Aerobee 350 sounding rocket, figure 5-1, are given in reference 14. 
Figures 5-2 and 5-3 show the results of the present method as well as data points from 
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Figure 5-1-Aerobee 350 





reference 14. The maximum derivation of the theoretical C N values is 7.83% and the maximum 
deviation of the theoretical center of pressure values is 1.18%. The degree of accuracy of the 
windtunnel tests is indicated as about 6%. 

5.20 ROLL FORCING MOMENT COEFFICIENT DERIVATIVE 

The only roll forcing moment coefficient derivative data over a broad range of mach numbers that are 
available to the author are those of the single stage tomahawk sounding rocket, figure 5-4. The data for this 
vehicle are unpublished at this time, but have been provided by the test coordinator, reference 28. 
Figure 5-5 shows the comparison between the theoretical result and windtunnel data. No data are 
available in the subsonic region. The average deviation of the theoretical C, is 7.54%. The 
accuracy of the windtunnel data is about 5%. It is apparent, both from the windtunnel data and the 
subsonic theory that the theoretical value at M = 1.5 is no good. 


5.30 ROLL DAMPING DERIVATIVE 

Reference 21 presents a full set of data for the basic finner test rocket, figure 5-6. It is the 
only configuration for which roll damping moment coefficient derivative data are available to the 
authors. Unfortunately, the range of the data is relatively small. Figure 5-7 shows the comparison 
between the theoretical c. is 5.68%. Reference 21 does not indicate the degree of accuracy of its 

P 

data. No subsonic data can be found for Q * 

p 

5.40 DRAG COEFFICIENT 

Aerobee 150A, figure 5-8, drag coefficient data are available from reference 29. The theoretical 
drag coefficient results are compared to the windtunnel data in figure 5-9. In order to duplicate 
the test conditions, the test reynolds numbers were used to set up a Re vs. M curve, figure 5-10. 

This curve was used in generating the theoretical data. The surface of the windtunnel model is 
highly finished. From table 4-1, the appropriate surface roughness is about .01 mils. The equiva- 
lent full scale surface roughness is .1 mils. The maximum deviation of the theoretical drag co- 
efficient is 8.52%. Again no subsonic data are available. 

5.50 PITCH DAMPING DERIVATIVE 

There are no c experimental data available. The data taken in windtunnel experiments are a 

mq 

combination of the pitch damping derivative and a time log derivative due to the rate of change of 
angle-of-attack. There is no way to separate the two effects. Therefore, no comparison of the 
pitch damping derivative can be attempted. 
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Figure 5-8-Aerobee 150A 
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Figure 5-9 — Aerobee 150A Drag Coefficient 
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6.00 CONCLUSIONS 


On the basis of the comparisons in section 5, the method presented in this dissertation pre- 
sented in this dissertation predicts the values of C Nt , X, , and to well within ten (10) 

percent of experimental results between mach 2 and 8. Although no comparison can be made for 
C mq , the results of section 5-1, coupled with the dependence of C mq on (C N , ) T(B) and X T(B) , tend to 
indicate that the theoretical c is also accurate to about 10%. 

mq 

The lack of subsonic experimental data prevents any concrete conclusions to be drawn about the 
accuracy of the present method at subsonic spfeeds. However, the figures in section 5 show that the 
subsonic variations of all the compared characteristics are quite reasonably valued with respect to 
the supersonic values. 

In general, then, the method presented in this dissertation can be used to confidently predict 
the aerodynamic characteristics of slender finned vehicles. 
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LIST OF SYMBOLS 


A Reference area 

C Length of airfoil strip chord 


Wave drag coefficient 


qA 


Lift coefficient - _ L 
qA 


dC. / da ( a = 0) 


Rolling moment coefficient = 


dC,/db ( = 0) 


qAL 


Pitching moment coefficient = 


Mp 
q AL 


C dC /da (3 = 0) 

m m 

a 

C Local pressure coefficient 

1 

C Length of root chord 

r 

d Component of 1 parallel to the freestream 

1 1 

d Portion of d behind the fin-tip Mach cone 

U W 1 

p Drag force (parallel to freestream) 

d 

F Lift force (normal to freestream) 

1 

K Interference factor 

L Reference length 

1 Length of airfoil-strip leading-edge region chord 

1 Length of root airfoil leading-edge region chord 

r 

1 Length of center region and of center-region chord 

1 Length of airfoil-strip trailing-edge region chord 

1 Length of root airfoil trailing-edge region chord 

r 

1 Distance between the leading edge and the intersection of the Mach cone with the airfoil 

strip 

1 Length of local region surface 

M Freestream Mach number 

M Pitching moment about the reference axis 
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LIST OF SYMBOLS (Continued) 



Rolling moment about the root chord 
Number of strips 

Component of L normal to the freestream 
Local static pressure 
Ambient static pressure 

% 

Freestream dynamic pressure 
d w /d ratio 
Semispan length 

Distance from the reference axis to the airfoil-strip leading edge in a streamwise direction 
Center-of-pressure coordinate measured from the reference axis 

Distance from the reference axis to the forward region boundary of a local region in a 
streamwise direction 

Center-of-pressure coordinate measured from the root chord 
Distance from the root chord in a spanwise direction 
Width of airfoil strip 
Angle of attack (degree) 

. \l 2 - 1 

Leading-edge sweep angle 
Trailing-edge sweep angle 
Ragion boundaries sweep angles 
Ratio of specific heats for air =1.4 
Fin cant angle 

Leading-edge wedge half-angle in the streamwise direction 
Trailing-edge wedge half-angle in the streamwise direction 
Inclination of a local surface to the freestream 
Mach cone semivertex angle 


SUBSCRIPTS 
B Body alone 

B(T) Body in the presence of the tail 

i Local region number 

j Airfoil strip number 
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LIST OF SYMBOLS (Continued) 


SUBSCRIPTS 
o Value at a = 0 

T Tail alone 

T(B) Tail in the presence of the body 
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FIN 


A COMPUTER PROGRAM FOR CALCULATING THE AERODYNAMIC 
CHARACTERISTICS OF FINS AT SUPERSONIC SPEEDS 


BACKGROUND 

Calculating the aerodynamic characteristics of the fins of a sounding rocket is one of the basic 
steps in calculating the overall aerodynamics of a sounding rocket. Because sounding rockets fly 
at supersonic speeds for the greater portion of their flight, the regime at these speeds is of particu- 
lar interest. 

Methods previously used for calculating supersonic fin aerodynamics (references 1, 2, 3, and 4) 
are lacking in accuracy and applicability because of the many assumptions and approximations in- 
herent in using them to reach a closed form or nearly closed form solution. 

By using a high-speed computer to numerically solve basic theoretical equations, one may ob- 
tain answers rapidly and as accurately as desired. The only restrictions to accuracy are the 
assumptions inherent in deriving the basic equations. The theory found to be most amenable to pro- 
gramming in Busemann's Second Order Supersonic Airfoil Theory as described in reference 1. 

PROGRAM CAPABILITIES AND LIMITATIONS 

Given a supersonic fin with chord plane symmetry, at a given Mach number, FIN computes as 
functions of angle of attack ( ) the following: 

• Lift coefficient, C L 

• Wave drag coefficient, C DW 

• Pitching moment coefficient, C m 

• Center-of- pressure coordinate measured from the reference axis, X 

• Center-of-pressure coordinate measured from the root chord, Y 

Since ; as defined in the program is equivalent to the fin cant angle ( : ) at i = 0, the rolling 
moment coefficient is considered a function of 

FIN also computes as functions of Mach number the following: 

• Lift coefficient slope, C L 

• Wave drag at zero angle of attack, C DW 

o 

• Pitching moment coefficient slope, C m 

m a 

• Rolling moment coefficient slope, 
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• Center-of-pressure coordinate measured from the reference axis at zero angle of attack, X 0 

• Center-of-pressure coordinate measured from the root chord at zero angle of attack, Y 0 

The range of is from zero to * max , as determined by the user, in increments of 1 degree. The 
program allows for a maximum of 100 Mach number points, the user specifying the range and 
increment of the Mach number. 

Since FIN has been designed for use with a four-finned vehicle, the output values of C L , C DB , 
r r c and C are for a pair of identical fins; the output values of C t and C J? are for two 

m ‘ J f m a » _ — . . , . j e. 

pairs of perpendicular fins. The values of X, Y, X Q , and Y Q apply only to a single typical fin. 

Busemann's Second Order Airfoil Theory is applied subject to the restrictions outlined in 
reference 1, pages 192 to 241, with special attention to Figures 10.3 and 10.28. However, the use 
of the third-order terms in this program slightly enlarges the applicability shown in this reference. 
Busemann's theory is applied to small streamwise airfoil strips, and the strips are summed in a 
spanwise direction. In addition, the fin-tip Mach cone is accounted for by applying a correction 
factor to the necessary portion of each strip. This correction-factor technique was obtained from 
references 1 and 3. The fin-tip Mach cone correction may or may not be included at the user’s 
discretion. 

Figure 1 shows the types of fins to which FIN is applicable. These configurations cover all 
of the shapes normally used on sounding rockets and missiles. The program input pattern for each 
type is given in the section on usage. A listing of FIN in FORTRAN IV is given in Appendix A. 

AERODYNAMIC THEORY 

Busemann’s Second Order Airfoil Theory has been applied to two-dimensional airfoil strips 
with the following assumptions: 

• All parts of the airfoil surface are in supersonic flow and make small angles with the flow. 
This implies low angles of attack. 

• The leading edge is sharp. The trailing edge must be sharp for a double wedge fin or a 
modified double wedge. For a single wedge and a modified single wedge, the effect of the 
blunt base is neglected. 

• The shock waves are all attached. 

• Each region of flow over the surface acts independently of the others. 

The basic local pressure coefficient equation of the Busemann theory is a third-order series 
expansion: 


C 


p 


P - P 

r i r °° 


K, Vi f *2^ + *3^ " 


( 1 ) 
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(o) Modified Double Wedgi 


(b) Double Wedgi 







3 


(c) Modified Single Wedge 


(d) Single Wedge 


Figure 1 — Fin Configurations 
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where 



If the flow over the surface region inclined at 77 . to the flow has been preceded by a single 
compressive shock wave anywhere upstream in the flow, the value of K* is: 



If the flow upstream of the region in question has contained only expansive shock-free flow, 

K* = 0. t, l is the inclination of the leading edge surface. Essentially, the first three terms of 
equation 1 represent the first-, second-, and third-order pressures on the surface region respec- 
tively. The K’ term is the third-order irreversible pressure rise through the bow shock. The 
index, i , indicates the surface for which the pressure coefficient is being calculated. 

The pressure coefficient of equation 1 is applied to a typical airfoil region, including the fin-tip 
Mach cone, by dividing the region into two portions at the intersection of the airfoil strip and the 
Mach cone (Figure 2). A correction factor of 1/2 is applied to the local pressure coefficient for the 
region within the Mach cone. The center of pressure of each portion is taken as the midpoint of 
the length of that portion. By this means, equations are obtained for the local lift ( F, _ ), drag ( F d . ), 
local center of pressure (x. ), and, hence, the local pitching moment (M p _ ): 
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Figure 2-Mach Cone Correction Geometry 
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The airfoil-strip lift (f, ), drag (F d ), and pitching moment (Mp ) totals are then computed by summing 
the local characteristics. 


V 1 

L 


(8) 



(9) 



(10) 


The strip rolling moment about the root chord is calculated from the airfoil lift and strip span- 
wise position, y : 

M r = yF, (11) 

j j 

The total characteristics over the entire fin are then computed by summing the airfoil strip 
characteristics in a spanwise direction. 


F 





F 


d 



M 


p 



M 

r 



( 12 ) 

(13) 

(14) 

(15) 


The forces and moments computed in equations 4 through 15 are actually divided by the ambient 
dynamic pressure (q^) and strip width (Ay). Thus, the associated coefficients are: 
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2F, • A .y 


( 2 fins) 


(16) 


2F - A.y 


"DW A 


(2 fins) 


(17) 


2M -Ay 

P J 

AL 


(2 fins) 


(18) 


4M r * Ay 
AL 


(4 fins) 


(19) 


The center-of-pressure coordinates are calculated from the respective moments and the 
appropriate lift forces. 


Fj cos 


( 20 ) 


( 21 ) 


To obtain the linear slope coefficients near a = 0 at a given Mach number, the values of C L , 
C m> and c, are calculated at a = 1° and are taken as the slope values c^, , and C, s (all per de- 

gree). The slope values per radian are obtained by the proper conversion of radians to degrees. 

Thus: 


C L 

J. 

rad. 

180 

77 

W. . ,) 

(22) 

C 

m 

a 

rad. 

180 
7 T 

('.1. . ,) 

(23) 

S 

r ad. 

180 

T 

( c .l. . ,) 

(24) 


The zero angle-of-attack values of C DW , X, and Y are simply their respective values at a - 0. 

CONFIGURATION AND GEOMETRY CONSIDERATIONS 

The user indicates the shape of the fin configuration by specifying the root chord length (C r ), 
root airfoil leading and trailing edge region chord lengths (l Lr and l Tr ), leading and trailing edge 


87 



sweep angles (I" L andr T ), region boundary sweep angles (r i and P 2 ), and the fin semispan length 
(s) (Figure 3). The program then computes the distance from the reference axis to the airfoil- strip 
leading edge in the streamwise direction (X L ), airfoil chord length (C), and the airfoil region chord 
lengths (l L and i T ) as functions of the spanwise coordinate of the airfoil strip. 


X L 

y 

t an i L cos x 


(25) 

c 

y 

( t an i T ~ t an • L ) 

* C 

r 

(26) 

>L 

z y 1 

( t an P , - t an f L ) 

+ i L 

r 

(27) 

>T 

= y! 

(tan P T - tan I ' 2 ) 

+ 1 T 

(28) 


The width of the airfoil strip is determined from the number of strips (n) desired by the user. 

.iy = — (29) 

n 

The airfoil strip crossection is assumed to be symmetrical about the chord line. The most 
general crossection configuration is the modified double wedge (Figure la). The other three types 
(Figures lb, lc, and Id) are special cases of this basic shape. Figure 4 shows in detail the general 
crossection, which has six separate flow regions. As can be seen in this figure, the surface region 
lengths are determined by specifying the leading- and trailing-edge wedge half angles in the stream- 
wise direction, £ and £ . 


(30) 


1 


3 


c - l L 


1 


T 


The program computes the local- surface inclination angles (r/. ) as functions of a: 

Tf i i ~ 7 
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(31) 

(32) 

(33a) 

(33b) 

(33c) 

(33d) 



SECTION A-A TYPICAL AIRFOIL 


Figure 3— Fin and Airfoil Geometry 





"7 


90 


'5 


(38e) 


% - ' T ♦ a (38f) 

From these, the streamwise and normal components of the region lengths are calculated. 


d j - 1 ( cos t ; | ( 1 1 - 6 ) 


n i 1 t ' s i n t { i : 1 - 6 ) 


(34) 


The distance of the forward local-region boundaries from the reference axis are then determined 
as: 


x p 1 

X FM X L 

(35a) 

X P 2 

X L * d , 

(35b) 

X P 3 

= X L * d, F d 2 

(35c) 

X P5 

X L ^ d < 

(35d) 

X P 6 

= X L ‘ <>« * d 5 

(35e) 


FIN-TIP MACH CONE CORRECTION 

The Mach-cone angle {&) is a direct function of the Mach number. 


Arc tan 



Ar cs in 



(36) 


For an airfoil strip at a spanwise direction y, the intersection of the strip and the Mach cone is a 
distance i w from the leading edge. 


l w (S - y)tan: L * t an ( 90 - (37) 

But, by trigonometric identities and the use of equation 37, this may be reduced to: 

l w (S - y) (tan ! \ ♦ ) (38) 

If is greater than or equal to the chord length, there is, of course, no actual intersection and no 
correction is necessary. If i is such that 


C 


l^C- 


91 


then the ratio of the portion of 1 falling behind the Mach cone to 1 T is: 


Also, 



If i w is such that 


then the corresponding ratio for the middle regions is 


Also, 


and, 


If l is such that 



l 


0 


C “ I T ” 1 M ’ U - 0 


then the ratio for the leading edge region is 


( 39 ) 


(40) 


(41) 


(42a) 

(42b) 


r 


1 



(43) 


Also, 


r 3 r 5 r 6 , -° 


(44) 


These ratios are the values used in equations 4, 5, and 6 for calculating the local lift, drag, and 
center of pressure. In these equations, the pressure coefficient in the portion behind the Mach cone 
is taken as half the value of the C p calculated by equation 1. 
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